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ABSTRACT 

The properties of the ices that form in dense molecular clouds represent an 
important set of initial conditions in the evolution of interstellar and preplane- 
tary matter in regions of active star formation. Of the various spectral features 
available for study, the bending mode of solid CO2 near 15 /xm has proven to be 
a particularly sensitive probe of physical conditions, especially temperature. We 
present new observations of this absorption feature in the spectrum of Q21-1, 
a background field star located behind a dark filament in the Cocoon Nebula 
(IC 5146). We show the profile of the feature be consistent with a two-component 
(polar + nonpolar) model for the ices, based on spectra of laboratory analogs with 
temperatures in the range 10-20 K. The polar component accounts for ~ 85% of 
the CO2 in the line of sight. We compare for the first time 15 /im profiles in three 
widely separated dark clouds (Taurus, Serpens and IC 5146), and show that they 
are indistinguishable to within observational scatter. Systematic differences in 
the observed CO2/H2O ratio in the three clouds have little or no effect on the 
15 /im profile. The abundance of elemental oxygen in the ices appears to be a 
unifying factor, displaying consistent behavior in the three clouds. We conclude 
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that the ice formation process is robust and uniformly efficient, notwithstand- 
ing compositional variations arising from differences in how the O is distributed 
between the primary species (H 2 0, C0 2 and CO) in the ices. 

Subject headings: Dust, extinction — Infrared: ISM: lines and bands 
— ISM: molecules 



1. Introduction 

Dust grains acquire icy mantles in the cold, dense molecular phase of the (ISM), the 
composition of which is determined both by gas-phase chemistry and by chemical reactions 
occurring on the surfaces of the grains (e.g., Tielens & Hagen 1982). Interstellar ices in- 
corporated into the envelopes and disks of young stellar objects (YSOs) may subsequently 
be heated, chemically or structurally modified or destroyed by energetic processes such as 
exposure to radiation, shocks and winds. The study of ices in the environments of YSOs is 
currently a highly active field, driven by the rich legacy of infrared spectra available from the 
Infrared Space Observatory (ISO) and the Spitzer Space Telescope (e.g., de Graauw et al. 
1996; d'Hendecourt et al. 1996; Gerakines et al. 1999; Nummelin et al. 2001; Gibb et al. 
2004; Boogert et al. 2004, 2008; Watson et al. 2004; Dartois et al. 2005; Pontoppidan et al. 
2005, 2008; Zasowski et al. 2008). Of the various infrared absorption features identified with 
vibrational modes of H 2 0, CO, C0 2 and other molecular species in the ices, the 0=C=0 
bending mode of C0 2 near 15 /im has proven to be particularly useful as a tool sensitive to 
physical conditions, especially temperature. Ices formed and maintained at low temperature 
are amorphous in structure, but studies of laboratory analogs have demonstrated the ap- 
pearance of subfeatures in the 15 /xm profile upon heating that are diagnostic of segregation 
and crystallization of constituents (Ehrenfreund et al. 1997, 1999; Dartois et al. 1999; van 
Broekhuizen et al. 2006; White et al. 2009). These subfeatures have, indeed, been observed 
in the infrared spectra of a variety of high- and low-mass YSOs (see Pontoppidan et al. 2008 
for numerous examples). 

The properties of pristine ices in cold, quiescent regions of dense molecular clouds pro- 
vide an important benchmark for studies of the evolution of interstellar matter and the 
origins of stars and planetary systems. To study ices in such regions, it is necessary to 
observe background field stars viewed through the dense material, and this is much more 
difficult compared with YSOs because of their relatively weak intrinsic flux at mid-infrared 
wavelengths. First detection of C0 2 -ice in a field star resulted from ISO observations of 
its near infrared (4.2 /im) stretching-mode feature in the spectrum of the Taurus K-type 
giant Elias 16 (Whittet et al. 1998). Subsequent observations with Spitzer confirmed the 
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presence of the b ending-mode feature in this line of sight (Bergin et al. 2005; Knez et al. 
2005), exhibiting a profile that is, as expected, notably lacking in structure associated with 
heating and crystallization of the ices. The feature is consistent with models that utilize data 
for laboratory analogs at low temperature (T ~ 10-20 K), combining 'polar' (H 2 0-rich) and 
'non-polar' (H^O-poor) components (e.g. Gerakines et al. 1999). Detections of CO2 in the 
spectra of a further 8 field stars behind the Taurus dark-cloud complex have since been 
reported (Whittet et al. 2007 and references therein). Consistency was found both in the 
CO2 abundance (~ 18% relative to H 2 0) and in the 15 /im profile shape: each line of sight 
displays a profile essentially identical to that of Elias 16 to within observational error. This 
result implies that the ices in Taurus have been maintained at low temperature since their 
formation, and furthermore that the formation process is fairly robust and uniformly efficient 
at several locations in the cloud. 

Taurus has (by default) become the prototype for study of ices in low-mass star- 
formation regions because of the relative ease with which they are detected in this dark-cloud 
complex0. A key question is whether it is a reliable prototype. Substantial cloud-to-cloud 
variations are known to exist, for example, in ambient radiation field (Liseau et al. 1999) and 
grain temperature (White et al. 1995), and also in gas phase chemistry between clouds of 
similar temperature, density and radiative environment (Dickens et al. 2000). It is important 
to determine whether such variations have an appreciable impact on the composition and 
physical state of the ices in order to test the validity of Taurus as a prototype. However, 
background field stars bright enough at 15 /im to yield spectra of sufficient quality and with 
enough extinction to show deep interstellar features are rare. To date the profile of the 
sensitive CO2 bending mode has been reported in only one confirmed background field star 
outside of the Taurus region: CK2 in Serpens (Knez et al. 2005) . The Serpens dark cloud dif- 
fers from that in Taurus in both structure and content: it is more centrally condensed, with 
higher average gas densities, and displays greater intermediate-mass star formation activity 
leading to warmer ambient gas temperatures (~ 25K compared with ~ 15K; White et al. 
1995; McMullin et al. 2000; Kaas et al. 2004). There is evidence also to suggest differences 
in ice composition between the two clouds, with (e.g.) an excess of CO and/or CO2 relative 
to H 2 in Serpens compared with Taurus (Eiroa & Hodapp 1989; Chiar et al. 1994; Whittet 
et al. 2007). 

In this paper, we report observations and profile modeling of the CO2 bending-mode ice 



1 The Taurus complex, which includes TMC-1, L1495, B18 and several other dense cores, is relatively 
nearby, with little foreground extinction. Ice features are detected at relatively low extinction compared 
with other clouds studied to date and are widely distributed across the region (see Shenoy et al. 2008 for 
further discussion and references). 
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feature in a field star background to a dark, filamentary cloud associated with the Cocoon 
Nebula (IC 5146) in Cygnus (the "Northern Streamer" mapped by Lada et al. 1994). We 
compare results for the three regions now available — Taurus, Serpens and IC 5146 - 
and show that they do, indeed, support the existence of a common 15 /im ice profile for 
unprocessed ices in the cold dark-cloud environment. Our results also allow new constraints 
to be placed on the ice formation process. 



2. Spectra 

All CO2 spectra used in this work were obtained with the Infrared Spectrograph (IRS) 
of the Spitzer Space Telescope, described by Houck et al. (2004). The IRS operated in 
short-wavelength, high-resolution staring mode, covering the spectral range 10.0-19.5 /xm at 
a resolving power A/ A A ~ 600. The Astronomical Observation Request (AOR) keys and 
other relevant information are given in Table 1 for our three program stars. The spectrum 
of Q21-1 (IC 5146), reported here for the first time, was obtained on 2004 November 12 as 
part of the Spitzer General Observer program #3320 (Principal Investigator: Y. Pendleton); 
three ramp cycles were carried out to yield a total on-source integration time of 7975 s, and 
a spectrum of the sky background in the region was also obtained to provide background 
subtraction. The raw spectra were processed by the standard Spitzer pipeline to produce 
basic calibrated data (BCD). Version 15.3.0 of the pipeline was used to reduce the data for 
Q21-1 and also for CK2 (a reduction of the CK2 spectrum using pipeline version 11.0.2 
was previously published by Knez et al. 2005). For Elias 16, we use the spectrum reported 
by Bergin et al. (2005), processed with version 13.0 of the pipeline^. Further analysis of 
the BCD to produce final calibrated spectra of Q21-1 and CK2 was carried out using the 
SMART program (Higdon et al. 2004) and standard procedures to remove bad pixels and 
other artifacts, as described in §2 of Whittet et al. (2007). 

The final calibrated flux spectrum of Q21-1 is shown in Fig. la. A local continuum was 
determined by fitting a second-order polynomial to adjacent wavelength regions on each side 
of the CO2 feature (13.0-14.9/im and 16.1-17.6/xm). The adopted continuum, shown as a 
dashed curve in Fig. la, was used to calculate the optical depth spectrum plotted in Fig. lb. 
The same procedure was used to produce an optical depth spectrum for CK2 (see below); 
our result for CK2 is broadly consistent with that reported by Knez et al. (2005), allowing 
for differences in the adopted continuum. 



2 Tests showed 110 improvement in data quality using pipeline version 15.3 to process the spectrum for 
Elias 16 compared with the published spectrum. 



- 5 - 



3. Results 

The final optical depth profiles of the CO2 features observed in CK2 and Q21-1 are 
shown in Fig. 2, together with fits based on laboratory analogs (Table 2). Fig. 3 (lower 
frame) compares the profiles of CK2 and Q21-1 with that of Elias 16, superposed with 
normalization to unity at peak optical depth. It is evident from inspection that the observed 
profiles in the three lines of sight are closely similar, indeed virtually identical to within 
observational error. This conclusion is supported by the random nature of the scatter in the 
residual optical depths for CK2 and Q21-1 relative to Elias 16, shown in the upper frame of 
Fig. 3. 

3.1. Comparison with laboratory analogs 

The fitting procedure used to model the CO2 profile is the same as described in our 
previous work (e.g. Gerakines et al. 1999; Bergin et al. 2005), but with access to a much 
larger suite of analog spectra now available. The procedure assumes that the profile can be 
modeled by a combination of two ice components^: one contains CO2 in a matrix dominated 
by H 2 (the polar component), the other contains CO and C0 2 , but little or no H 2 (the 
nonpolar component). Some analogs also contain methanol (CH3OH), which can affect 
the shape of the C0 2 bending mode (Ehrenfreund et al. 1999). A total of 1381 laboratory 
spectra is available from the databases of the Astrophysics Laboratories of Leiden University 
(Ehrenfreund et al. 1997, 1999) and the University of Alabama at Birmingham (White et al. 
2009), spanning a wide range of temperature (5-150 K) as well as composition. 

We fit each spectrum independently to further test the hypothesis that the profiles 
are effectively the same. Best fits resulting from our least-squares fitting procedure are 
summarized in Table 2 and plotted in Fig. 2 for CK2 and Q21-1 (our fit to Elias 16 is 
similar to that shown in Fig. 2 of Bergin et al. 2005). As expected, the combinations of 
ices selected are qualitatively and quantitatively similar, with temperatures in the range 10- 
20 K; indeed the same polar ice mixture was selected in each case, and the nonpolar mixtures 
differ only in CO:C0 2 ratio. In fact, the fits are not very sensitive to C0 2 concentration in 
either component — different mixtures give almost equally good fits (see also §4.1 below) - 
but they are robust in terms of three essential requirements: (i) both polar and nonpolar 
components must be present, (ii) both components must be at low temperature (< 25 K), 



3 The choice of two as the number of components is justified as the smallest capable of giving an acceptable 
fit (no single mixture can reproduce the profile), but it should be noted that as many as five components 
have sometimes been used to fit CO2 bending-mode profiles in YSOs (Pontoppidan et al. 2008). 
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and (iii) the profile is dominated by the polar component, which contributes ~ 80-90% of 
the total absorption in each line of sight. None of the CH 3 OH-bearing mixtures were selected 
by the fitting process, consistent with available limits on its solid-state abundance in dark 
clouds (Chiar et al. 1996). 

3.2. Column densities 

The integrated optical depth, J r(u)du, of the observed absorption feature in each pro- 
gram star was estimated, and results are listed in Table 1. This quantity is directly related 
to column density: 

N = t^ (1) 
A 

where A is the band strength. Gerakines et al. (1995) showed that the band strength of the 
C0 2 bending mode is only weakly dependent on the ice mixture in which the C0 2 resides, 
and we therefore adopt an average value of A = 1.0 x 10~ 17 cm/molecule, consistent with our 
previous work (Whittet et al. 2007). Column densities for C0 2 are listed in Table 2, together 
with values for other major constituents of the ices (H 2 and CO) from the literature. All 
iV(H 2 0) values are determined from ground-based observations of the 3.0 /zm O-H stretching 
mode, which generally gives more reliable results than the blended O-H-0 bending mode 
near 6 fim. 

Fig. 4 plots CO and C0 2 ice column densities against that of H 2 0. In addition to 
Elias 16, CK2 and Q21-1, all Taurus field stars with available data are included (Whittet 
et al. 2007 and references therein). The straight lines are linear least-squares fits to Taurus 
stars only. The fit to iV(CO) vs. iV(H 2 0) displays a positive intercept on the iV(H 2 0)- 
axis, as expected, as H 2 is less volatile and can exist in solid form at lower dust columns 
compared with CO. It is evident (Fig. 4a) that both CK2 in Serpens and Q21-1 in IC 5146 
have significantly higher C0 2 /H 2 abundance ratios compared with members of the Taurus 
group (which show a rather tight correlation). Systematic differences are less apparent in the 
plot of iV(CO) vs. iV(n 2 0), suggestive of a common trend (Fig. 4b) albeit with appreciable 
scatter; the locus of Q21-1, in particular, is closely consistent with the Taurus correlation. 

4. Discussion 

The results described in the previous section demonstrate that the C0 2 ice feature near 
15 /iin displays a common profile shape in three widely separated dense clouds. The available 
spectra are consistent with a combination of polar (H 2 0-rich) and nonpolar (CO-rich) ices 
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at low temperature (10-20 K), the polar component accounting for ~ 85-90% of the CO2. 
Any elevation in ice temperature arising from proximity to embedded stars (e.g., Serpens: 
see §1) must not be sufficient to affect the profile in the lines of sight studied. Interestingly, 
this commonality of profile also occurs notwithstanding a clear difference in the composition 
of the ices, manifested by discontinuity in the iV(C02) vs. iV(H20) correlation for IC 5146 
and Serpens compared to Taurus (Fig. 4a). In this section, we explore possible implications 
of these results. We hypothesize that the compositional discontinuity between clouds may 
arise from variations in only one parameter: the CO2/H2O ratio in the polar ices. 

4.1. Stability of the profile 

It was noted previously that the bending-mode absorption profiles of ices containing 
CO2 are not strongly sensitive to CO2 concentration at a given temperature. This remains 
true provided only that the ices maintain their compositional and structural 'identity' (e.g., 
a polar mixture contains enough H 2 to remain 'polar'). For illustration, we compare in 
Fig. 5 the profiles for several polar laboratory ice mixtures at low temperature, with CO2 
concentration varying between 14% and 100% relative to H2O. Slight differences in position 
and width are seen, but the profile shape is maintained over this rather broad compositional 
range. It is therefore unsurprising that observed profiles are closely similar in lines of sight 
with different CO2/H2O ratios. 

4.2. Oxygen budget 

Polar ices are thought to form on interstellar grains by surface reactions involving species 
(principally H, O and CO) adsorbed from the gas phase in the "diffuse-dense" transition 
region of a molecular cloud where H — > H 2 conversion is incomplete. In brief, the main 
chemical pathways appear to be hydrogenation of O to OH and H2O, and oxidation of CO 
to CO2, with no more than trace amounts of the CO undergoing hydrogenation to form 
CHO-bearing species such as CH3OH (e.g., Tielens & Whittet 1997)5 The C0 2 /H 2 ratio 
in the resulting solid is likely to depend on the relative abundances of the adsorbed species, 
and thus on the network of chemical reactions occurring in the gas phase. Local conditions 
may lead to variations. A cloud with a relatively high concentration of gas phase CO relative 



4 Note that the chemistry leading to CO2 formation is not yet fully understood because of uncertainties in 
the rates of key surface reactions such as CO + O — > CO2 and CO + OH — > CO2 + H; see Ruffle & Herbst 
(2001) and Pontoppidan et al. (2008) for further discussion. 
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to atomic H in the transition region will tend to produce a more C0 2 -rich polar ice, whereas 
a cloud with a lower concentration of gas phase CO will produce less C0 2 and free up more 
O to form H 2 0: the differences evident in Fig. 4a between Taurus and the other clouds might 
be explained in this way. At higher densities, H — > H 2 conversion is essentially complete and 
H 2 can no longer form efficiently by surface reactions; nonpolar ices then accumulate on the 
grains, principally by attachment of gas-phase CO (Pontoppidan 2005), which is abundant 
in the dense molecular phases of all dark clouds. Formation of C0 2 within the nonpolar 
component may be driven by cosmic-ray processing of the solid CO (Pontoppidan et al. 
2008). 

It is informative to consider the distribution of elemental oxygen in the ices arising from 
these processes. Ices, together with silicate dust and gas-phase CO, represent the major 
observable oxygen reservoirs in the dense ISM (Whittet et al. 2007 and references therein). 
O is the most abundant heavy element and the only one present in all three of the major 
molecular ice species, H 2 0, CO and C0 2 . Atomic O attaching to a grain surface will undergo 
chemical reactions that lead to formation of H 2 and C0 2 , as outlined above: we hypothesize 
that the fraction of O atoms accumulating into the ices is robust but that their subsequent 
distribution between H 2 and C0 2 may vary from cloud to cloud. 

The total number of O atoms in the polar ices along a line of sight may be expressed in 
terms of column density by the relation 

iV(0) polar = iV(H 2 0) + hN (CO) + 2/ 2 iV(C0 2 ), (2) 

where f\ and / 2 are the fractions of all solid-phase CO and C0 2 molecules, respectively, in 
the polar ices. Similarly, for the nonpolar ices, we have 

iV(0) nonpolar = (1 - A)JV(CO) + 2(1 - / 2 )JV(C0 2 ). (3) 

Previous and current work suggests f\ ps 0.17 in Serpens and 0.12 in Taurus and IC 5146 
(Chiar et al. 1994, 1995 and in preparation); / 2 is evaluated from the results of the present 
paper (Table 2, right hand column). We ignore other species that might contain significant 
elemental oxygen in the ices, of which the most important is probably 2 (see Whittet et al. 
2007 for further discussion). 

Fig. 6 plots N(0) for polar and nonpolar ices against visual extinction (Ay). As in 
Fig. 4, all available Taurus field stars are included using data from Whittet et al. (2007) in 
addition to results from the present work (Elias 16 has the highest A v value of the Taurus 
stars plotted). Ay is taken as a measure of total dust column. High degrees of correlation 
are evident in Fig. 6, with notable consistency between the three clouds for both polar 
and nonpolar ices: in contrast to the situation in Fig. 4, a single linear fit to all points is 
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representative of the trend for each component of the ices. This result suggests stability in 
the net rate of accumulation of O atoms into ices on grain surfaces, regardless of differences 
in the efficiencies of chemical pathways that form H 2 and C0 2 . 

Some caveats should be noted. First, the extinction toward CK2 is particularly un- 
certain: we have adopted the value from Knez et al. (2005), but higher values have been 
estimated in earlier literature (see discussion in §6 of Whittet et al. 2007). Second, we 
are concerned only with the component of Ay arising in dense clouds (where grains are 
ice-mantled); any systematic cloud-to-cloud differences in the extinction contributed by un- 
mantled grains in diffuse phases of the ISM along these lines of sight would lead to systematic 
error in the form of a horizontal shift in Fig. 6. This contribution may be estimated empir- 
ically by evaluating the ice 'threshold' extinction, i.e., the intercept on the Ay axis of fits 
to correlations of H 2 0-ice optical depth or column density with visual extinction: a value of 
~ 3.3 mag is well documented for the Taurus cloud (Whittet et al. 1988, 2001), and this is, 
indeed, consistent with the locus of the intercept of the linear fit to the polar component in 
Fig. 6. Available data suggest somewhat higher thresholds for IC 5146 (~ 4.0 mag; Chiar 
et al. in preparation) and Serpens (~ 5.5 mag; Eiroa & Hodapp 1989). However, in each 
case, the difference relative to the Taurus value is smaller than the estimated uncertainty in 
Ay for Q21-1 and CK2, respectively. Therefore, we do not consider this potential source of 
error to be of major significance. 

Fig. 6 suggests a uniform distribution of elemental oxygen between polar and nonpolar 
components, the latter accounting for ~ 20% of the total in the ices at high Ay. This 
conclusion is supported by Fig. 7, which plots iV(O) for polar ices against that for nonpolar 
ices (eqs. 2 and 3). Fig. 7 is unaffected by any systematic error introduced by cloud-to- 
cloud difference in diffuse-ISM extinction, and may be compared with Fig. 4b (the largest 
contributions to O in polar and nonpolar ices being made by H 2 and CO, respectively). The 
data are consistent with a linear correlation with continuity between clouds. The positive 
intercept of the fit on the iV(0) po iar axis in Fig. 7 is probably real (c.f. Fig. 4b) and consistent 
with the lower volatility of polar ices compared with nonpolar ices. 

4.3. Homogeneity of polar/nonpolar ice phases 

Consistency in the distribution of matter between polar and nonpolar phases of the 
ice is indicated both by the correlations in Figs. 6 and 7 and by the results of the 15 /zm 
profile fitting (§3.1). An unexpectedly high level of compositional homogeneity is implied: 
no major, systematic variations in the relative abundance of polar to nonpolar ice is apparent 
over a rather wide range of environments and extinctions (5 < Ay < 35 mag; Fig. 6). Polar 
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ices form at relatively low density, as discussed above, and their abundance is expected to 
correlate linearly with dust column (i.e., with Ay) if mantled grains are widely distributed 
throughout the clouds. In the case of nonpolar ices, however, the abundance of the primary 
species (CO) is not dependent on formation in situ but rather on accumulation from the gas. 
The rate of accumulation would thus be expected to increase rapidly with density at high 
Ay unless there is an efficient desorption mechanism. Our results suggest that a balance 
is maintained between accretion and desorption of CO. This situation is, indeed, entirely 
consistent with observations of gas-phase CO in molecular clouds in general (substantial 
levels of gas phase CO are maintained in dense regions), and with laboratory work on the 
photodesorption properties of CO ice (see Oberg et al. 2007 for further discussion). 



5. Conclusions 

The primary result of this study is a clear demonstration that the C0 2 ice bending- 
mode feature near 15 /im displays essentially the same profile shape in lines of sight that 
intercept three widely separated dense clouds. This profile is likely to be a useful benchmark 
for investigations of C02-bearing ices in YSOs, aiding separation of contributions made by 
the circumstellar environment and the parent molecular cloud to each observed spectrum. 

The shape of the observed profile elucidates the nature of the ices in the dense-cloud 
environment. Comparisons with laboratory analogs confirm that the C0 2 -bearing ices are 
maintained at low temperature (< 25 K) and include both polar (H 2 0-dominated) and non- 
polar (CO-dominated) components. The polar component contains ~ 85-90% of all solid- 
phase C0 2 , but the concentration of C0 2 relative to H 2 is not tightly constrained. The 
constancy of the profile in the face of observed differences in the C0 2 /H 2 abundance ratio 
is thus easily understood. 

Cloud to cloud variations in C0 2 /H 2 ratio may result from differences in gas phase 
chemistry, but the physical process of element depletion from the gas onto dust proceeds 
in a predictable manner. The elemental oxygen column density in the ices appears to be 
a unifying parameter, showing evidence of consistent behavior with respect to dust column 
in our sample. Moreover, the distribution of elemental oxygen between polar and nonpolar 
ice appears to be almost invariant: the CO-rich nonpolar component accounts for no more 
than ~ 20%, even in the densest regions sampled, a result which implies the existence of an 
efficient desorption mechanism for CO in molecular cloud cores. 

It will be important in the future to extend this work to a larger sample of field stars 
that probe quiescent regions of molecular clouds. Given the constancy of profile observed to 



- 11 - 



date, spectra obtained with the low-resolution mode of the Spitzer IRS may be adequate to 
extend abundance correlations such as those shown in Figs. 4 and 6, allowing fainter stars 
to be included. Opportunities to obtain high-quality spectra of other ice constituents best 
observed from airborne and/or space platforms will be provided in the next few years by the 
availability of the Stratospheric Observatory for Infrared Astronomy and the James Webb 
Space Telescope. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under a con- 
tract with NASA. The authors acknowledge NASA support from grant NNX07AK38G and 
JPL/Caltech Support Agreement no. 1290823 (D.C.B.W.) and JPL/Caltech subcontracts 
1266411 and 1267778 (J.E.C.). We are grateful to an anonymous referee for helpful com- 
ments that led to improvements to this paper. 
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Table 1. Summary of observations and related information for program stars. 



2MASS ID 


Other ID 


Cloud 


Spectral 


A v 


AOR key 




Notes b 








type 


(mag) 




(cm -1 ) 




J04393886+2611266 


Elias 16 


Taurus 


Kl III 


24.1 ± 0.5 


03868160 


5.4 ±0.2 


(1) 


J18300061+0115201 


CK2, EC 118 


Serpens 


K4 III 


34 ±4 


11828224 


9.9 ±0.4 


(2) 


J21472204+4734410 


Q21-1 


IC 5146 


K2 III 


27 ±1 


10751488 


8.7 ±0.4 


(3) 



integration was carried out over the wavelength range 14.9-16.1 £tm. 

b Notcs: (1) IRS spectrum from Bergin et al. 2005; spectral type from Elias 1978; visual extinction (Ay) from 
Shenoy et al. 2008. (2) IRS spectrum previously reported by Knez et al. 2005 and rereduced by us (see §2); 
spectral type and extinction from Knez et al. 2005. (3) IRS spectrum from the current work (§2); spectral type 
and extinction from Chiar et al., in preparation. 



Table 2. Column densities and details of fits to the CO2 bending-mode feature. 



Star 


Column densities" 


Component 


Mixture 


T 


Percent 




iV(H 2 0) 


iV(CO) 


JV(C0 2 ) 






(K) 


N(C0 2 ) 


Elias 16 


25.7 b 


6.5 C 


5.4 


polar 


H 2 0:C0 2 :CO = 100:20:3 


20 


89% 










nonpolar 


CO:C0 2 = 100:23 


10 


11% 


CK2 


35 d 


16.2° 


9.9 


polar 


H 2 0:C0 2 :CO = 100:20:3 


20 


88% 










nonpolar 


CO:C0 2 = 100:26 


10 


12% 


Q21-1 


25.2 f 


7.1 f 


8.7 


polar 


H 2 0:C0 2 :CO = 100:20:3 


20 


85% 










nonpolar 


CO:CQ 2 = 100:70 


10 


15% 



a All column densities are in units of 10 17 cm 2 ; all Af(C0 2 ) values are from the current work. 
b Average value from 3 fim spectra of Whittet et al. 1988, Smith et al. 1993 and Murakawa ct al. 2000. 
c Chiar ct al. 1995. 

d Calculated from 3/j.m spectrum of Eiroa & Hodapp 1989. 

e Average of values from Chiar et al. 1994 and Pontoppidan ct al. 2003 (as cited by Knez et al. 2005). 
'From unpublished observations (Chiar et al., in preparation). 
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Fig. 1. — (a) Flux spectrum of Q21-1 (IC 5146) in the spectral range 13.0-17.6 fim. The 
continuum adopted for the CO2 feature centered near 15.2 jum is shown as a dashed curve, 
(b) Optical depth spectrum over the same spectral range, obtained using the continuum 
shown in (a). 
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Fig. 2. — Optical depth spectra of (a) CK2 and (b) Q21-1, with our best-fitting model to the 
C0 2 profile superposed in each case: polar component (dashed curves), nonpolar component 
(dotted curves), total (thick solid curves). The mixtures used in the fits are listed in Table 2. 
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Fig. 3. — Below: superposition of observed C0 2 profiles in all three field stars, normalized 
to unity at peak optical depth (Elias 16: thick solid curve; CK2: thin dotted curve; Q21-1: 
thin solid curve). Above: residuals obtained by subtracting the normalized profile of Elias 16 
from those of CK2 and Q21-1 (dotted and solid curves, respectively). 
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Fig. 4. — Correlations of column densities for ices in field stars: (a) A^(H 2 0) vs. iV(C02) and 
(b) iV(H20) vs. iV(CO). Solid circles: Elias 16 and other Taurus stars from Whittet et al. 
(2007); open squares: Q21-1 (IC 5146); open diamonds: CK2 (Serpens). The diagonal line 
is the linear least-squares fit to Taurus stars only in each case. The slopes and intercepts 
(respectively) of the fits are (a) 0.179 ±0.042, 0.055 ±0.554; (b) 0.342 ±0.026, -1.558 ±0.450. 
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Fig. 5. — Superposition of normalized CO2 bending-mode profiles in laboratory ices at low 
temperature with differing H 2 0:C02 ratios. The thick curve represents the polar mixture 
used in the fits shown in Fig. 2 (H 2 0:C0 2 :CO = f00:20:3 at 20 K); the other curves all 
represent H 2 0:C0 2 mixtures at fO K: 1:1 (dot-dash), 5:1 (dashed), and 100:14 (dotted). 
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Fig. 6. — Plots of the estimated column density N(0) of elemental oxygen in polar and 
nonpolar ice components (filled and open symbols, respectively) against the total visual 
extinction Ay. Taurus, Serpens and IC 5146 stars are represented by circles, diamonds and 
squares, respectively. The diagonal lines are linear least-squares fits to the polar and nonpolar 
data sets, combining results for all three clouds. The slopes and intercepts (respectively) 
of the fits are 0.175 ± 0.035, -6.646 ± 0.625 (polar), and -2.883 ± 0.764, 0.485 ± 0.043 
(nonpolar). 
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Fig. 7. — Plot of the estimated column density N(0) of elemental oxygen in polar ice vs. 
that in nonpolar ice. Symbols have the same meaning as in Fig. 4. The diagonal line is the 
linear least-squares fit to the combined data set for all three clouds. The slope and intercept 
of the fit have values 0.293 ± 0.036 and -1.774 ± 1.090, respectively. 



